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Abstract 

Lie symmetry group method is applied to study for the higher order Camassa-Holm equation. The symmetry group and its 
optimal system are given. Furthermore, preliminary classification of its group invariant solutions, symmetry reduction and 
nonclassical symmeries are investigated. Finally conservation laws for the higher order Camassa-Holm equation are presented. 
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j^. | 1 Introduction 

t-H ■ 

^ In the study of shallow water waves, Camassa and Holm [7] derived a nonlinear dispersive shallow water wave 
equation 

Ut — u x 2 t + 3uu x = 2u x u x 2 + uu x 3, (1) 

which is called Camassa-Holm equation (CH). Here u(x,t) denotes the fluid velocity at time t in the x direction or, 
k>( \ equivalently, the high of water's free surface above a fiat bottom. Eq. (1) has a bi-Hamiltonian structure [13,17], and 
■ is completely integrable [7,11]. It has many conservation laws [21]. Moreover, the CH equation is a re-expression of 
' geodesic flow on the diffeomorphism group of the line. Holm, Marsden and Ratiu [14] have shown that CH equation 
in n dimensions describes geodesic motion on the diffeomorphism group of R n with respect to metric given by the 
H 1 norm or Euclidean fluid velocity. Misiolek [19] has shown that the CH equation represents a geodesic flow on the 
Bott-Virasoro group. Kouranbaeva [15] has shown that the CH equation (for the case k = 0) is a geodesic spray of 
the weak Riemannian metric on the diffeomorphism group of the line or the circle obtained by the right translation 
of the H 1 inner product over the entire group. This equation admits well known properties and a rich literature is 
devoted to it. 



In recent years, many researchers have been researched on the Camassa-Holm equation. They extend the studies 
to the generalized CH equation, higher order CH equations and so on. Lixin Tian, Chunyu Shen and Danping 
Ding gave the optimal control of the viscous CH equation under the boundary condition and proved the existence 
and uniqueness of optimal solution to the viscous CH equation in a short interval (see [26]). Using geometrical 
methods, higher order CH equations have been treated in [12]. The well-posedness of higher order CH equations 
were considered in [10]. 



Email addresses: m_nadjafikhaii@iust.ac.ir (Mehdi Nadjafikhah), v.shirvani@kiau.ac.ir (Vahid Shirvani-Sh) . 



Preprint submitted to Automatica 



23 November 2011 



The formulation of the higher order Camassa-Holm equation which was recently derived by Coclite, Holden and 
Karlsen in [10] is 



B k (u,u) := A^Ckiu) -uu. 



A k (u) := £(-l)^j U , 



(2) 



C k (u) = -uA k {d x u) + A k (ud x u) - 2d x uA k (u). 

where k is a positive integer. In cases k = and k = 1, Eq. (2) becomes the inviscid Burgers equation and the 
Camassa-Holm equation respectively. 

In this paper we only consider the case k = 2 of Eq. (2). It also can be rewritten as 



The theory of Lie symmetry groups of differential equations was developed by Sophus Lie [16], which was called 
classical Lie method. Nowadays, application of Lie transformations group theory for constructing the solutions 
of nonlinear partial differential equations (PDEs) can be regarded as one of the most active fields of research in 
the theory of nonlinear PDEs and applications. Such Lie groups are invertible point transformations of both the 
dependent and independent variables of the differential equations. The symmetry group methods provide an ultimate 
arsenal for analysis of differential equations and is of great importance to understand and to construct solutions 
of differential equations. Several applications of Lie groups in the theory of differential equations were discussed in 
the literature, the most important ones are: reduction of order of ordinary differential equations, construction of 
invariant solutions, mapping solutions to other solutions and the detection of linearizing transformations. For many 
other applications of Lie symmetries see [22,6,4]. 

The fact that symmetry reductions for many PDEs arc unobtainable by applying the classical symmetry method, 
motivated the creation of several generalizations of the classical Lie group method for symmetry reductions. The 
nonclassical symmetry method of reduction was devised originally by Bluman and Cole in 1969 [5], to find new exact 
solutions of the heat equation. The description of the method is presented in [8,18]. Many authors have used the 
nonclassical method to solve PDEs. In [9] Clarkson and Mansfield have proposed an algorithm for calculating the 
determining equations associated to the nonclassical method. A new procedure for finding nonclassical symmetries 
has been proposed by Bila and Niesen in [1]. 

Many PDEs in the applied sciences and engineering are continuity equations which express conservation of mass, 
momentum, energy, or electric charge. Such equations occur in, e.g., fluid mechanics, particle and quantum physics, 
plasma physics, elasticity, gas dynamics, electromagnctism, magneto-hydro-dynamics, nonlinear optics, etc. In the 
study of PDEs, conservation laws are important for investigating integrability and linearization mappings and for 
establishing existence and uniqueness of solutions. They are also used in the analysis of stability and global behavior 
of solutions [2,3,24,25]. 

This work is organized as follows. In section 2 we recall some results needed to construct Lie point symmetries of a 
given system of differential equations. In section 3, we give the general form of a infinitesimal generator admitted by 
equation (3) and find transformed solutions. Section 4, is devoted to the nonclassical symmetries of the higher order 
CH model, symmetries generated when a supplementary condition, the invariance surface condition, is imposed. 
In Section 5, we construct the optimal system of one-dimensional subalgebras. Lie invariants, similarity reduced 
equations and differential invariants corresponding to the infinitesimal symmetries of equation (3) are obtained in 
section 6. Finally in last section, the conservation laws of the equation (3) are obtained. 

2 Method of Lie Symmetries 

In this section, we recall the general procedure for determining symmetries for any system of partial differential 
equations see [22,20,6,4]. To begin, let us consider the general case of a nonlinear system E of partial differential 
equations of order n in p independent and q dependent variables is given as a system of equations 



A := u t - u x 2 t + u x i t + 3uu x - 2u x u x 2 - 



UU X 3 + 2u x U x 4 + UU x a = 0. 



(3) 



A„(xy")) = 0, v = !,-■■, I, 



(4) 
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involving x — (x 1 , • • • , x p ), u — (it 1 , • • • , u q ) and the derivatives of u with respect to x up to n, where represents 
all the derivatives of u of all orders from to n. We consider a one-parameter Lie group of infinitesimal transformations 
acting on the independent and dependent variables of the system (4) 

x i =x i + sC(x,u) + 0{s 2 ), i = l---,p, ^ 

v? = it- 7 + s(p J (x, u) + 0(s 2 ), j = 1 • ■ • , q, 

where s is the parameter of the transformation and if are the infinitesimals of the transformations for the 
independent and dependent variables, respectively. The infinitesimal generator v associated with the above group of 
transformations can be written as 

P Q 

v = Y,C{x,u)d xi + Y,V> a {^u)d u c (6) 

i— 1 a—1 

A symmetry of a differential equation is a transformation which maps solutions of the equation to other solutions. 
The invariance of the system (4) under the infinitesimal transformations leads to the invariance conditions (Theorem 
2.36 of [22]) 

Pr (n) v[A„(a:,u( n >)] = 0, v = l,---,l, whenever A„(x, u (n) ) = 0, (7) 
where 

Pr («) 

is called the n th order prolongation of the infinitesimal generator given by 
Pr(")v = v+^^^( a; ,nW)5 u «, (8) 

a=l J 

where J = (ji, ■ ■ ■ ,jk), 1 < jk < P, 1 < k < n and the sum is over all J's of order < #J < n. If #J = k, the 
coefficient ip" of d u ^ will only depend on fc-th and lower order derivatives of u, and 

ipiix, = Dj(<p a - co + E (9) 
»=i »=i 

where uf := du a /dx l and Uj 4 := duj/dx 1 . 

One of the most important properties of these infinitesimal symmetries is that they form a Lie algebra under the 
usual Lie bracket. 

3 Lie symmetries for the higher order CH equation 

We consider the one parameter Lie group of infinitesimal transformations on (x 1 = x, x 2 = t, u 1 = u), 

x = x + s£(x, t, u) + 0{s 2 ), 

t = x + sr](x,t,u) + 0(s 2 ), (10) 
U = X + SLf(x, t, u) + 0(s ), 

where s is the group parameter and £ x — £, £ 2 = r\ and ip 1 = if are the infinitesimals of the transformations for the 
independent and dependent variables, respectively. The associated vector field is of the form: 

v = £(x, t, u)d x + r](x, t, u)d t + ip(x, t, u)d u . (11) 

and, by (8) its fifth prolongation is 

Pr^v = v + <p x d Ux + p l d Ut + p x2 d Ux2 + p xt a Uxt + p t2 d Ut2 

+p* 3 d Ux3 + p* 2t d Ux2t + p xt2 d Uxt2 + p t3 d Ut3 +■■■ + p xt4 d Uxt4 + d Ut5 . (12) 
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where, for instance by (9) we have 

ip x = D x (<p - £u x - t]u t ) + £u x 2 + rju xt , 
= D t (ip - £ u x - 77 u t ) + S, u xt + 77 u t 2 , 

; (13) 

</ = D l{^-i u x-'nut)+£,u xH +riu t $, 

where D x and D t are the total derivatives with respect to x and t respectively. By (7) the vector field v generates a 
one parameter symmetry group of the Eq. (3) if and only if 

Pr( 5 V[ut — u x 2 t + u x t t + 3uu x — 2u x u x 2 — uu x 3 + 2u x u x 4 + uu x s] = 0, whenever A = 0. (14) 
The condition (14) is equivalent to 

(3u x - u x 3 + u x s)tp + ip t + (3u - 2u x 2 + 2u x i)ip x - 2u x ip x2 - wp* 3 - Lp^ 1 + 2u x ip xi + uip x ° + Lp x ^ 1 = 0, 

(15) 

whenever w t — u x 2 t + u x <t t + 3uu x — 2u x u x 2 — uu x 3 + 2u x u x i + uu x s = 0. 

Substituting (13) into (15), and equating the coefficients of the various monomials in partial derivatives with respect 
to x and various power of u, we can find the determining equations for the symmetry group of the Eq. (3). Solving 
this equations, we get the following forms of the coefficient functions 

£ = c 3 , ri = cit + c 2 , (p = -ciu. (16) 

where c\, C2 and C3 are arbitrary constant. Thus, the Lie algebra g of infinitesimal symmetry of the Eq. (3) is spanned 
bye the three vector fields 

vi=<9 x , v 2 =<9 t , v 3 =td t -ud u . (17) 

The commutation relations between these vector fields are given in the Table 1. The Lie algebra g is solvable, because 
if g (1) = (vj, [vj, vj]) = [0,0], we have (1) = (vi,v 2 ,v 3 ), and = [0 (1) ,0 (1) ] = (v 2 ), so, we have a chain of ideals 

Table 1 

The commutator table 



[Vi,Vj] 


Vl 


v 2 


V3 


VI 











v 2 








v 2 


V3 





-V 2 






To obtain the group transformation which is generated by the infinitesimal generators v, for i = 1, 2, 3 we need to 
solve the three systems of first order ordinary differential equations 



ds 
dt(s) 

ds 
du(s) 



■■7]i(x(s),t(s),u(s)), t(0)=t, 1-1,2,3 
■■ipi(x(s),t(s),u(s)), u(Q) = u. 



ds 

Exponentiating the infinitesimal symmetries of equation (3), we get the one- parameter groups Gi(s) generated by 
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v, for i = 1, 2, 3 



G 2 
G 3 



(t,x,u) i 

(t, X, u) ! 

(i,a;,ti) i 



(x + s,t,u), 
(x, t + s, u), 
(x, e%e~ s u). 



(18) 



Consequently, 

Theorem 3.1 If u = f(x,t) is a solution of higher order CH equation, so are the functions 

G 1 (s)-f(x,t) = f(x-s,t), 
G 2 (s)-f(x,t) = f(x,t-s), 
G 3 (s)-f(x,t) = f(x,te- s )e- s . 



(19) 



Nonclassical symmetries for the higher order CH equation 



In this section we would like to apply the nonclassical method to the higher order CH equation. The graph of a 
solution 

u a = f a (x u ---,x p ), a = l,---,q (20) 

to the system (4) defines an p-dimensional submanifold Ff C R p x H q of the space of independent and dependent 
variables. The solution will be invariant under the one-parameter subgroup generated by vector (6) if and only if 
Tf is an invariant submanifold of this group. By applying the well known criterion of invariance of a submanifold 
under a vector field we get that (20) is invariant under vector (6) if and only if / satisfies the first order system Eq 
of partial differential equations 

p 

Q a (x,u,u^)=^ a \x,u)-J2C^^>i=0, a = 1, •••,<? (21) 

i=l 

known as the invariant surface conditions. The q-tuple Q = (Q 1 , • • • , Q q ) is known as the characteristic of the vector 
field (6). In what follows, the n-th prolongation of the invariant surface conditions (21) will be denoted by Eq\ 
which is a n-th order system of partial differential equations obtained by appending to (21) its partial derivatives 
with respect to the independent variables of orders j < n — 1. 

For the system (4), (21) to be compatible, the n-th prolongation Pr^v of the vector field v must be tangent to the 

(n) 

intersection E ("1 E k q ' 

Pr(")v(A,)| £nB („ ) =0, i/ = l, (22) 

If the equations (22) are satisfied, then the vector field (6) is called a nonclassical infinitesimal symmetry of the 
system (4). The relations (22) are generalizations of the relations (7) for the vector fields of the infinitesimal classical 
symmetries. A similar procedure is applicable to the case of the nonclassical infinitesimal symmetries with an 
evident difference that in general one has fewer determining equations than in the classical case. Therefore, we 
expect that nonclassical symmetries are much more numerous than classical ones, since any classical symmetry is 
clearly a nonclassical one. The important feature of determining equations for nonclassical symmetries is that they 
are nonlinear, this implies that the space of nonclassical symmetries does not, in general, form a vector space. For 
more theoretical background see [23,1]. 

If we assume that the coefficient of d t of the vector field (6) does not identically equal zero, then for the vector field 
v = t, u)d x + d t + (p(x, t, u)d u (23) 



5 



the invariant surface conditions are 



ut + £,u x = ip, (24) 

Calculating equations (22) and inserting ip from (24) in to it, we can find the determining equations by equating the 
coefficients of the various monomials in partial derivatives with respect to x and various power of u. Solving this 
equations, we get £ = c and f — 0, where c is arbitrary constant. 

Now assume that the coefficient of <9 t in (23) equals zero and try to find the infinitesimal nonclassical symmetries of 
the form 

v = d x + tp(x, t, u)d u (25) 

for which the invariant surface conditions is u x = tp. Similar the previous case, we can find determining equations. 
Solving this equations, we get ip = 0. This means that no supplementary symmetries, of non-classical type, are 
specific for our models. 



5 Optimal system for the higher order CH equation 

In general, to each s-parameter subgroup H of the full symmetry group G of a system of differential equations in 
p > s independent variables, there will correspond a family of group-invariant solutions. Since there are almost 
always an infinite number of such subgroups, it is not usually feasible to list all possible group-invariant solutions to 
the system. We need an effective, systematic means of classifying these solutions, leading to an "optimal system" of 
group-invariant solutions from which every other such solution can be derived. 

Definition 5.1 Let G be a Lie group with Lie algebra q. An optimal system of s— parameter subgroups is a list of 
conjugacy inequivalent s— parameter subalgebras with the property that any other subgroup is conjugate to precisely 
one subgroup in the list. Similarly, a list of s— parameter subalgebras forms an optimal system if every s— parameter 
subalgebra of g is equivalent to a unique member of the list under some element of the adjoint representation: 
f, = A%(f,)), fl eG.[22] 

Theorem 5.2 Let_H and H be connected s- dimensional Lie subgroups of the Lie group G with corresponding Lie 
subalgebras f) and f) of the Lie algebra g of G. Then H=gHg~ 1 are conjugate subgroups if and only if f) = Ad(g(i))) 
are conjugate subalgebras. (Proposition 3.7 of [22]) 

By theorem (5.2), the problem of finding an optimal system of subgroups is equivalent to that of finding an optimal 
system of subalgebras. For one-dimensional subalgebras, this classification problem is essentially the same as the 
problem of classifying the orbits of the adjoint representation, since each one-dimensional subalgebra is determined 
by nonzero vector in g. This problem is attacked by the naive approach of taking a general element V in g and 
subjecting it to various adjoint transformation so as to "simplify" it as much as possible. Thus we will deal with the 
construction of the optimal system of subalgebras of g. To compute the adjoint representation, we use the Lie series 

Ad(exp(evi)vj) = vj - e[vi,\j] + —[vi, [v i7 vj]] , (26) 

where [vj, Vj] is the commutator for the Lie algebra, e is a parameter, and i,j = 1, 2,3. Then we have the Table 2. 
Table 2 



Adjoint representation table 



Ad(exp(evi)vj) 


Vl 


v 2 


V3 


Vl 


Vl 


v 2 


V3 


V2 


Vl 


V2 


V3 — £V2 


V3 


Vl 


V2 + £V 3 


V3 



Theorem 5.3 An optimal system of one- dimensional Lie algebras of the higher order CH equation is provided by 
(1) avx + vz, (2) fivi + v 2 
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proof: Consider the symmetry algebra g of the Eq. (3) whose adjoint representation was determined in table 2 and 
let F? : q — > q defined by Ad(exp(eVj)v) is a linear map, for i = 1,2, 3. The matrices Mf of Ff , i = 1,2, 3, with 
respect to basis {vi,v 2 ,v 3 } are 



M{ 



/ioo\ 

1 

V° V 



M 9 £ 



(l o o\ 

1 e 
1 



MS 



/l o\ 

e~ £ 

V° i/ 



Let V = 5Z i=1 a i v « i s a nonzero vector field in q. We will simplify as many of the coefficients ai as possible by acting 
these matrices on a vector field V alternatively. 

Suppose first that a 3 ^ 0, scaling V if necessary we can assume that a 3 — 1, then we can make the coefficients of v 2 
vanish by Mf, and V reduced to case 1. 

If a 3 = and a 2 ^ 0, then we can not make vanish the coefficients of vi and v 2 by acting any matrices MS. Scaling 
V if necessary, we can assume that a 2 = 1 and V reduced to case 2. □ 



6 Symmetry reduction and differential invariants for the higher order CH equation 



Lie-group method is applicable to both linear and non-linear partial differential equations, which leads to similarity 
variables that may be used to reduce the number of independent variables in partial differential equations. By 
determining the transformation group under which a given partial differential equation is invariant, we can obtain 
information about the invariants and symmetries of that equation. 

Symmetry group method will be applied to the (3) to be connected directly to some order differential equations. 
To do this, a particular linear combinations of infinitesimals are considered and their corresponding invariants 
are determined. The equation (3) is expressed in the coordinates (x, t, u), so to reduce this equation is to search 
for its form in specific coordinates. Those coordinates will be constructed by searching for independent invariants 
(y, v) corresponding to the infinitesimal generator. So using the chain rule, the expression of the equation in the 
new coordinate allows us to the reduced equation. Here we will obtain some invariant solutions with respect to 
symmetries. First we obtain the similarity variables for each term of the Lie algebra g, then we use this method to 
reduced the PDE and find the invariant solutions. 



We can now compute the invariants associated with the symmetry operators, they can be obtained by integrating 
the characteristic equations. For example for the operator aw\ + v 2 = ad x + d t characteristic equation is 

dx dt du 

~ = T = T' ( ' 

The corresponding invariants are y = x — at, v = u therefore, a solution of our equation in this case is u = v{y). 
The derivatives of u are given in terms of y and v as 

Ut = —av y , u x 2 t = -av y 3, u x 4 t = -av y 5, u x — v y , u x 2 = v y 2, ■ ■ ■ , u x s = v y a. (28) 

Substituting (28) into the Eq. (3), we obtain the ordinary differential equation 

— av y + a v y 3 — a v y s + 3 v v y — 2 v y v y 2 — v v y 3 + 2 v y v y * + v v y s = 0. (29) 

All results are coming in the tables 3 and 4. 

Differential invariants help us to find general systems of differential equations which admit a prescribed symmetry 
group. One say, if G is a symmetry group for a system of PDEs with functionally differential invariants, then, the 
system can be rewritten in terms of differential invariants. For finding the differential invariants of the equation (3) 
up to order 2, we should solve the following systems of PDEs: 

di di di di 

W (30) 
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Table 3 



Reduction of Eq. (3) 


operator 


y 


V 


u 


Vl 


t 


u 


v{y) 


v 2 


X 


u 


v(y) 


V3 


X 


t u 




a vi + V3 


x - log(i) 


tu 


Hv) 


OV1+V2 


x — at 


u 


v{y) 



Table 4 

Reduced equations corresponding to infinitesimal symmetries 



operator 



similarity reduced equations 



3VVy 



2VyV y 2 



vv y 3 + 2v y v y 4 + VVyS 







Vl 

v 2 

V3 — V + Vy2 — Vyi + 3VVy — 2VyVy2 — VVy3 4" 2V y V y 4 4" VVy5 = 

a vi + V2 —a Wj, + a w y 3 — av y 5 +3vv y — 2 v y v y 2 — v v y 3 +2 v y v y i +vv y 5 =0 

a vi + V3 —v — av y + v y 2 + a v y 3 — v y i — a v y s + 3vv y — 2v y v y 2 — vv y 3 + 2v y v y i + vv y s = 



where / is a smooth function of (x, t, u), 

dh dh 



8h 

dx ' 



dh 

dt ' 



dh_ 

•x ^ 

ou x 



dt du 

where I\ is a smooth function of (x, t, u, u x ,u t ), 



dh 

dx ' 



dh 
dt ' 



t 



dh 
dt 



dh 
du 



dh_ 

du x 



2m, 



2u t 



dh 

du t ' 



dh 
du t 



dh 
1x2 « — 

du T : 



2u xt 



dh 

du x 



3 u t 2 



dh 

duv. 



(31) 



(32) 



where h is a smooth function of (x 7 t,u 7 u x ,Ut 7 u XX7 u x t,u t t)- The solutions of PDEs systems (30), (31) and (32) 
coming in table 5, where * and ** are refer to ordinary invariants and first order differential invariants respectively. 

Table 5 

differential invariants 



vector field 


ordinary invariant 


1st order 


2nd order 


Vl 


t, u 


*,u x ,u t 


*, **, Uxx-, Uxt, Utt 


v 2 


X, u 


*,u x ,u t 


*, **, U xx , U xt , Utt 


v 3 


x,tu 


*,tu x ,t 2 u t 


*,**,tu xx ,t 2 U xt ,t 3 Utt 



7 Conservation laws for the higher order CH equation 

Many methods for dealing with the conservation laws are derived, such as the method based on the Noether's 
theorem, the multiplier method, by the relationship between the conserved vector of a PDE and the Lie-Bcklund 
symmetry generators of the PDE, the direct method, etc. [22,2,3,24]. 

Now, we derive the conservation laws from the multiplier method. 

Definition 8.1 A local conservation law of the PDE system (4) is a divergence expression 

A$>] = D^iu] + ■■■ + D r ^ n [u] = (33) 

holding for all solutions of the system (4). In (33), = <& l {x, u, d u , • • • , <9£), i = 1, • • ■ , n, are called fluxes of the 

conservation law, and the highest-order derivative (r) present in the fluxes [u] is called the order of a conservation 
law. [3] 
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In particular, a set of multipliers {A„[[/]}J, =1 = {A„(x, U, du, • • • , du)}l=i yields a divergence expression for the 
system A v (x, u'"') such that if the identity 

A V [U]A V [U] = Di&[U] (34) 

holds identically for arbitrary functions U(x). Then on the solutions U(x) = u(x) of the system (4), if A„[[/] is 
non-singular, one has local conservation law Ai,[u]A„[u] = = 0. 

Definition 8.2 The Euler operator with respect to CP is the operator defined by 

^=^7-A^ 7 + - + (-irA 1 -A 5 ^ + - (35) 

tl'"l s 

for j = 1,-.. ,q.[3] 

Theorem 8.3 The equations E u , F(x,U,du, ■ ■ ■ ,d v ) = 0, j — I,--- ,q hold for arbitrary U{x) if and only if 
F{x, U, du, • • • , dfj) = D i ^! t (x, U, du, • • • , d^ 1 ) holds for some functions ^ % {x, U, du, ■■■ , d^' 1 ), i = 1, • ■ • q. (The- 
orem 1.3.2, [3]) 

Theorem 8.4 A set of non-singular local multipliers {A v (x, U,du, ■ • ■ , d\j)} l v=1 yields a local conservation law for 
the system A v (x,u^ n ') if and only if the set of identities 

Eui (A„(x, U,du,--- , d r u)A„(x, U W)) = 0, j = 1, • • • q, (36) 
holds for arbitrary functions U(x). (Theorem 1.3.3, [3]) 

The set of equations (36) yields the set of linear determining equations to find all sets of local conservation law 
multipliers of the system (4). Now, we seek all local conservation law multipliers of the form A = £(x, t, u) of the 
equation (3). The determining equations (36) become 

Eu[Z(x,t, U)(U t - U xH + U xH + 3UU X - 2U X U X 2 - UU x s + 2U X U X * + UU X ,)] = 0, (37) 

where U (x, t) are arbitrary function. Equation (37) split with respect to third order derivatives of U to yield the 
determining PDE system whose solutions are the sets of local multipliers of all nontrivial local conservation laws of 
the higher order CH equation. 

The solution of the determining system (37) given by 

i = Cl U + c 2 , (38) 
where c\ and ci are arbitrary constants. So local multipliers given by 

1K = 1, 2)£ = U, (39) 

Each of the local multipliers £ determines a nontrivial local conservation law D t ^> + D x <fr = with the characteristic 
form 

A* + D x § ee i{U t - U xH + U xH + 3UU X - 2U X U X 2 - UU x s + 2U X U X * + UU x s), (40) 

To calculate the conserved quantities ^ and $, we need to invert the total divergence operator. This requires the 
integration (by parts) of an expression in multi-dimensions involving arbitrary functions and its derivatives, which is 
a difficult and cumbersome task. The homotopy operator [25] is a powerful algorithmic tool (explicit formula) that 
originates from homological algebra and variational bi-complexes. 
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Definition 8.5 The 2-dimcnsional homotopy operator is a vector operator with two components, (j^^ x t )/i ^( x t)f) ' 



dX 



(41) 



where 

4(V=/ 1 (E« , /)[A«]^ and til, t) f=^(±$f)[Xu] 
The x-integrand, / ^ /, is given by 

Mf M 2 3 fel _i fe 2 

#/= E E ( E E ^^(-d^-- 1 (-a)*-* )-^, (42) 

fe 1= lfc 2= i 1= 0i 2 =0 OU x k it k 2 

where , M| are the order of / in u to x and t respectively and combinatorial coefficient 



B (I) =B(ii, i 2 ,fci,fc 2 ) = 



/il + *2\ /«1 + fc 2 - «1 - 12 - 1\ 

\ ii / \ fci - h - 1 / 



Similarly, the t-integrand, /*• /, is defined as 

"<*.*) 

Mf m| fel fe 2 _i 

$/= E E (E E ^^^(-^^(-a)*--"- 1 )^-. (43) 

fe 1= 0fc 2 = l i!=0i 2 =0 OU x k lt k 2 

where B^(i 2 , h, k%, k\). 

We apply homotopy operator to find conserved quantities <J/ and $ which yield of multiplier £ = 1. We have 

f=U t - U xH + U xH + 3UU X - 2U X U X 2 - UU x3 + 2U X U X * + UU X , , (44) 
the integrands (42) and (43) are 

/ = 3 u 2 - u 2 x -\u xt -2 uu x 2 + 2 u x u x 3 - u 2 x2 + f u xH + 2 uu x i , 

(45) 

4 t * ) /= u-\u x 2 + \u x *, 
apply (41) to the integrands (45), therefore 

* := ^u(x,t)f = \u 2 -\u 2 x -\u xt - uu x 2 + u x u x 3 - \ u 2 x2 + I U x 3 t + UU x i, 

(46) 

so, we have the first conservation low of the higher order CH equation respect to multiplier £ = 1 

Ae(§ u 2 - \ u 2 x - | u xt - uu x 2 + u x u x 3 -\u 2 £ 2 + \ u x 3 t + uu x i) + D t (u - ± u x 2 + ± u^) = 0. (47) 
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Now we find conservation law respect to multiplier £ = u, in this case we have 

/= U(U t - U xH + U xH + 3UU X - 2U X U X 2 - UU x s + 2U X U X , + UU X ,), (48) 
the integrands are 

/f/= 3m 3 + 3m 2 v - \ uu xt + \u x u t -3u 2 u x 2 + § uu x s t - ^u x su t - \ u x u xH + \ u x iu xt , 

(49) 

/*■/= u 2 - | uu x 2 + §«£ + § uu x i - § u x u x s + \u 2 x2 , 
applying 2-dimensional homotopy operator, we have 

* := ^u(x,t)f = ^ + W 2 ^ 4 - I UU xt + I U x U t - U 2 U X 2 + \UU XH - \ U x 3U t - I U x U x 2 t + I U x 2U xt , 

(50) 

* := / 4*(x,t)/= 5 " 2 - I Wi ^ 2 + i " 2 + I uu ^ 4 - I u ^ w ^ 3 + To u x 2 ' 
so, the second conservation low of the higher order CH equation is 

D x (u 3 + U 2 U X 4 - | UU xt + | U x U t - U 2 U X 2 + ^UU x 3 t - | U^s^ - | U x U x 2 t + | U x 2U xt ) 

(51) 

" 2 - \ uu x 2 + \ u l + \ uu x* - \ u x u x 3 + ± u 2 x2 ) = 0. 
8 Conclusion 

In this paper by applying the criterion of invariance of the equation under the infinitesimal prolonged infinitesimal 
generators, we find the most general Lie point symmetries group for higher order CH equation. By applying the 
nonclassical symmetry method for the higher order CH equation, we concluded that the analyzed model do not 
admit supplementary, nonclassical type symmetries. Also, we have constructed the optimal system of one-dimensional 
subalgebras for the higher order CH equation. The latter, creates the preliminary classification of group invariant 
solutions. The Lie invariants and similarity reduced equations corresponding to infinitesimal symmetries are obtained. 
We find the conservation laws from the multiplier method. 
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